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ABSTRACT 

An analytical investigation of aeroelastic response of metallic and com- 
posite propfan models in yawed flow has been performed. The analytical model 
is based on the normal modes of a rotating blade and the three-dimensional 
unsteady lifting surface aerodynamic theory including blade mistuning. The 
calculated blade stresses or strains are compared with published wind tunnel 
data for two metallic and three composite propfan wind tunnel models. The com 
pari son shows a good agreement between theory and experiment. Additional para 
metric results indicate that blade response is very sensitive to the blade 
stiffness and also to blade frequency and mode shape mistuning. From these 
findings, it is concluded that both frequency and mode shape mistuning should 
be included in aeroelastic response analysis. Furthermore, both calculated 
and measured strains show that combined blade frequency and mode shape mistun- 
ing has beneficial effects on response due to yawed flow. 


INTRODUCTION 

Renewed interest in propfan (or advanced turboprop) propulsion is due to 
its potential for very high fuel efficiency at cruise speeds up to Mach 0.85. 

A number of studies have been conducted by both NASA and industry to evaluate 
the potential of propfan propulsion for both civil and military applications. 
These studies have led to two promising propfan concepts: One is single rota- 

tion (SR) with eight to ten blades and the other is counterrotation (CR). 

Both these concepts have been extensively studied. Large scale models have 
been built, and flight-tested (refs. 1 and 2) for proof of concept. 
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Supporting these new propfan concepts, both theoretical and experimental 
research effort on propfan aeroelasticity has been continuing at NASA Lewis 
Research Center. As a part of this research, a theoretical investigation on 

forced aeroelastic response of SR propfans in yawed- or off-axis flow has been 

performed. The yawed flow causes a periodic change in the angle of attack of 
blades, and this change, in turn, causes one per revolution periodic loading 
on the blades. Understanding and predicting this kind of periodic loading and 
the propfan blade response are very important in the design of the blades. 

Theoretical prediction of propfan aeroelastic response is one of the most 

challenging problems faced by aeroelasticians because the blades are thin and 
flexible with large sweep. The deflections due to centrifugal and aerodynamic 
loads are large, particularly at the blade tip; hence the response problem is 
inherently nonlinear requiring geometric nonlinear theory of elasticity in 
structural modeling of the blades. Since the blades have a moderate aspect 
ratio with large sweep, three-dimensional aerodynamic effects in addition to 
cascade effects should be included in the aerodynamic modeling of the blades. 
Additionally, the properties of the blades differ from blade to blade because 
of manufacturing limitations. The presence of such small differences is known 
as mistuning and it is expected to play a significant role on the blade aero- 
elastic response. Because of these complexities, published literature on the 
propfan forced response is very limited. References 3 and 4 calculated the 
forced response of three propfan wind tunnel models tested at NASA Lewis. 
Analysis was based on two-dimensional unsteady aerodynamic theory without 
blade mistuning. The correlation between theory and experiment varied from 
poor to good and identified a need for improved analysis methods. 

The specific objectives of the present investigation are: (1) to improve 

the propfan forced response prediction capabilities by developing a theoretical 
model, and an associated computer program, with three-dimensional unsteady 
aerodynamic effects and aerodynamic and structural mistuning; (2) to correlate 
theoretical and experimental results for validating the analysis; and (3) to 
conduct parametric studies to assess the effect of frequency and mode shape 
mistuning on blade response. 

To accomplish these objectives, the theoretical aeroelastic models 
described in references 5 and 6 were enhanced to calculate blade stresses due 
to yawed flow. As in references 5 and 6, each blade was structurally modeled 
by using plate finite elements. The three-dimensional, unsteady subsonic lift- 
ing surface aerodynamic theory, described in reference 7, was extended to 
account for blade mistuning and for yawed flow. The aerodynamic theory was 
also extended to cover supersonic helical Mach numbers with subsonic axial 
flow. By using blade normal modes and the associated unsteady aerodynamic 
loads, a modal forced response analysis and an associated computer program, 
ASTR0P3, were developed. The acronym ASTR0P3 stands for, "Aeroelastic Stabil- 
ity and Response of Propulsion Systems." 

The code was applied to five wind tunnel propfan research models of 0.62 m 
<2 ft) diameter which were tested for forced response at NASA Lewis. These 
models are SR-5 (ref. 8), SR-3 (ref. 4) SR3C-X2 (ref. 5), SR3C-3 (ref. 6), and 
SR3C-X2/SR3C-3 (ref. 6). The blades of the first two models were made of tita- 
nium, and of the last three models were molded from graphite-ply/epoxy-matrix 
laminated material. The blade mistuning in the first four models was inherent 
from manufacturing tolerances and was relatively small. However, the fifth 
model was intentionally mistuned for a particular type of mistuning which is 
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characterized by the differences in the blade frequencies, mode shapes, steady- 
state deflected shapes, and motion dependent unsteady aerodynamic loads. 

The paper presents the analytical development of an aeroelastic response 
model, outlines the application of the analysis to the wind tunnel models, and 
correlates calculated and measured blade stresses. Finally, the paper presents 
some results to show the effects of mistuning on the aeroelastic response of 
the propfan models. 
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NOMENCLATURE 

matrix of influence coefficients, equation (5) 
generalized aerodynamic matrix 
area vector, equation (6) 
matrices, equations (10) and (11) 
nonaerodynamic force vector of a group of blades 
modal frequency of a blade in jth mode, Hz 
generalized aerodynamic force vector due to yawed flow 
unit matrix 

also blade index 

generalized stiffness matrix of a group of blades 
generalized stiffness matrix ith blade in a group 
stiffness matrix of a group of blades, equation (1) 
centrifugal softening matrix of a group of blades 
aerodynamic matrix, r = 0,1 . . . NG - 1 , equation (8) 
Mach number of freestream 
physical mass matrix 

generalized mass matrix of a group of blades 
generalized mass matrix of a jth blade 
generalized mass of ith blade in jth mode 
number of blade groups 
number of blades in a group 


3 


NMG 

NP 

NMi 

NRT 

[P] 

{P({u 0 }>} 

{q} 

{q 0 > 

j 


[Si] 

{u 0 } 

[H] 

&0.75R 

{Au} 

[S] 

5iJ 

a Gi 

a xxi , a yyi >°xyi 

C*1] 

w 1j 


number of modes of a group 
number of aerodynamic panels in a group 
number of modes of the ith blade in a group 
total number of panel edge radii in a group 
aerodynamic loads of a blade group 

steady-state aerodynamic modal force vector of a group of 
blades 

generalized coordinate vector of a blade group 
complex amplitude of {q} 

generalized coordinate of the jth mode of the ith blade 
in a group 

modal stress matrix of an element of the ith blade 

steady-state deflection 

matrix of normal velocities 

blade pitch angle at 3/4 radius 

vibratory deflection vector at grid points 

normal displacement of the blade surface at the control points 

damping ratio in jth mode of the ith blade in a group 

intergroup phase angle between like blades from group to group 

normal and shear stresses of an element of the ith blade 

modal matrix of ith blade in a group 

frequency of the jth mode of the ith blade, rad/sec 


FORMULATION OF ANALYTICAL MODEL FOR AEROELASTIC RESPONSE 

An analytical model for flutter and forced response by the modal method 
for tuned single-rotation propfans with rigid hubs was presented in refer- 
ence 5. That model was extended and generalized in reference 6 for flutter of 
single rotation propfans with structurally and aerodynamical ly mistuned blades. 
Effects of the aerodynamic mistuning on the propfan flutter (including alter- 
nately mounted backward and forward swept blades) were addressed in refer- 
ence 9. In this section, the model presented in reference 6 has been extended 
and generalized for aeroelastic response of SR propfans in yawed flow. Since 
the basic details were given in references 5 and 6, only an outline of the 
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development and of special modifications which are essential to incorporate 
both aerodynamic and structural mistuning are presented here. 


Equations of Motion for Aeroelastic Response of a Mistuned Propfan 

The propfan, shown in figure 1, is assumed to have NG identical groups 
of blades symmetrically distributed around the disk. Each group contains NB 
blades which need have no special geometric or spatial relationship to each 
other. A schematic for an eight-bladed rotor with four blade groups is illus- 
trated in figure 2. A tuned rotor then corresponds to NB = 1, in which case 
NG is the total number of blades. In general, it is assumed that the groups 
vibrate with identical motion, but with a constant intergroup phase angle, 
which may take any of the values crgi = 2-iri/NG, i = 0,1, ... (NG-1). General- 
izing the steady-state single blade equation in reference 5, the corresponding 
equation for any group is 

[[K $ ] + [K<{u 0 }>]] {u Q } = {P({u 0 })} + {F o } (1) 

The form of equation (1) is the same as the corresponding one in refer- 
ence 5, but the orders of the matrices and of the displacement vector are dif- 
ferent. The stiffness and force matrices in equation (1) represent the entire 
blade group. For example, for a propfan with eight blades (four groups with 
two blades in each group), 228 grid points for each blade, and six degrees of 
freedom for each grid point, the total number of degrees of freedom of equa- 
tion (1) is 2736. Furthermore, the equation is nonlinear involving large 
deflections, centrifugal loads, and steady-state airloads. The solution 
requires substantial computer time. To reduce that time, the steady-state 
aerodynamic loads have been neglected herein (the effect of these steady air- 
loads on flutter for a tuned rotor is addressed in ref. 5). With this 
assumption, the steady aerodynamic coupling between the blades in a group is 
eliminated from the analysis. Consequently, the steady-state deflections and 
the differential stiffness of each blade in a group can be calculated indepen- 
dently by using equation (1) without steady airloads. The computational pro- 
ceure used for both the metallic and composite blades is the same as that 
described in reference 5. The vibration modes and frequencies of each blade 
in the group are calculated independently by using the following equation, 
which is also discussed in reference 5. 

[M] {Aii(t)} + [[K $ ] + tK({u 0 })]] {Au( {t) } = 0 (2) 

For calculating the blade elemental stresses and strains, the blade modal 
stresses are required. In solving equation (2), the modal stresses are also 
calculated. The next step in the formulation of the aeroelastic model by the 
modal method is to generalize the single blade equation of reference 5 to NB 
blades in a group. Assuming simple harmonic motion and following the procedure 
in reference 6, the aeroelastic equations for NB arbitrarily mistuned blades 
in a group are 

£ Mg g {q} + t K g g {q} = [A,] {q} + {G}e 1ut (3) 

Again, the form of equation (3) is the same as the corresponding one of 
reference 5, but the orders of the matrices and of the generalized coordinate 
vector are different. The new definitions are as follows: 
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The coefficients M^j, toi-j, and Cij are the generalized mass, the fre- 
quency, and the damping ratio, respectively, corresponding to the jth normal 
mode of the ith blade in a group. The suffix NMi represents the number of 
blade normal modes considered for the ith blade in the analysis. The quanti- 
ties C<t>i ] and (q i > are the modal matrix and the generalized coordinate 
vector for the ith blade in a group, respectively. The details for the 
development of the motion dependent aerodynamic matrix, [A]], and the motion 
independent aerodynamic matrix, {G}, are given in the following section. 

Unsteady Aerodynamic Model For Arbitrary Mistuning 

The generalized aerodynamic force matrix, [Ail, in equation (3) has been 
computed with the three-dimensional compressible lifting surface theory 
described in reference 7. Since the extension of the theory to include mistun- 
ing is given in reference 6, only an outline of the method and special modifi- 
cations for calculating the aerodynamic matrix, {G> , will be discussed. 

The numerical method is based on an assumed linearization of the fluid 
motion about a uniform steady flow. The disturbance generated by the rotor 
blades is presumed to be simple harmonic in time in a reference frame rotating 
with the blades. Each wake is assumed to be rigid and to lie in the helical 
surface swept out by the blade trailing edge. These assumptions result in a 
linear integral equation relating the normal velocity of the surface of the 
blades to the load (differential pressure) acting on the blades. The integral 
equation is solved by a piecewise constant load panel method. 

For simplicity, the effects of hubs and nacelles are ignored, as is blade 
thickness. The blades are, in effect, treated as loaded helical surfaces with 
no thickness. As in any linear aerodynamic theory, the unsteady loads induced 
by blade vibration are mathematically decoupled from the steady loads arising 
from blade incidence. 

The use of linear aerodynamic theory limits the applicability of the 
method to lightly loaded conditions (unstalled) and to Mach numbers at which 
transonic effects at the blade tips are not significant. In addition, the 
rigid wake approximation breaks down at low advance ratios, where the induced 
velocities may have a strong influence on the wake structure (which in turn may 
significantly modify the blade loading). Of course, one usually does not know 
apriori whether such effects will be present at a given operating point, so 
some caution in interpreting the results is necessary. 

The geometry and vibration mode data are usually obtained by solving 
equations (1) and (2) with COSMIC NASTRAN (ref. 10). For aerodynamic model- 
ing, the group is paneled by placing a fixed number, NXP, of quadrilateral pan- 
els on each radial interval r t ( j ) < r < rt(j + 1), j = 1, ... NRT (unless 
rt(j + 1) < rt(j), which signals a jump from the tip of one blade to the root 
of the next). Thus, except for the deletion of tip-root connections, the group 
is effectively treated as a single surface with NP = NXP* (NRT - NB) panels. 

The definition of group vibration modes requires special mention. Essen- 
tially, one group mode consists of one blade in the group vibrating in one of 
its natural modes while all other blades in the group remain fixed. Suppose 
the jth blade in a group is assigned NMj vibration modes. The group as a 
whole, then has NMG = NMI + NM2 + ... + NMNB modes. The group mode NMI + 1, 
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for example, corresponds to the second blade in the group vibrating in its 
first mode, while all remaining blades are stationary. 

In general, let matrix [W(n,k>] denote the normal velocity at the nth 
control point, the point at which the boundary condition is enforced, in the 
kth group mode. This array is initialized to zero and then loaded sequen- 
tially with the NM1 values of first blade modes [<t>i 3 , at control points on 
the first blade, the NM2 values of second blade modes, [> 2 ], at control 
points on the second blade, and so on. The normal velocities at each control 
point due to the yawed flow are stored in the last column of the [W] matrix. 
Thus the total number of rows of the [W] matrix is equal to the number of 
control points in the group, and the total number of columns is equal to 
NMG + 1 . 


The aerodynamic loads can then be found by solving the linear algebra 
problem: 


[A] [PJ = [H] 


(5) 


where [A] is the (NP,NP) matrix of influence coefficients, [W] is the (NP, 
NMG + 1) matrix of normal velocities, and [P3 is the (NP, NMG + 1) array of 
loads on the blade group. The influence coefficient matrix contains the 
interference between blade groups and, therefore, depends on the intergroup 
phase angle (as well as the vibratory frequency). Thus, separate solutions 
are required for every intergroup phase angle and frequency of interest. For 
example, for an eight-bladed propfan with tuned blades, the intergroup phase 
angle of the forcing function is 315°, which corresponds to a one nodal diam- 
eter backward traveling wave. For an eight-bladed propfan with four symmetri- 
cally identical groups, the intergroup phase angle of the loading is 270°. 

Once the load distribution has been found, the generalized force matrix 
can be determined by numerical integration over the blade group: 

NP 

A^n.m) = P(j,m) $(j,n) dA(j) (6) 

3=1 

where P(j,m) is the load on the jth panel (with area dA(j)), associated 
with the mth group mode, and $(j,n) is the normal displacement of the sur- 
face at the control point of the jth panel in the nth mode. Also, the 
aerodynamic force matrix due to yawed flow is 

NP 

G(n> = X P<;3, NMG + ]> 8<;3 ’ n) dA( ^ <7) 

j = l 

The matrices [Ai ] and {G} can be calculated for any given frequency 
of oscillation and any given intergroup phase angle. However, for yawed flow 
the frequency of the aerodynamic loading on the blade is the same as the rota- 
tional frequency of the propfan and the intergroup phase angle of the loading 
is 2ir(NG - 1)/NG. 


The above formulation (eq. (3)), is valid for an arbitrarily mistuned 
propfan. Also in this case, the entire rotor has to be modeled for calculating 
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the aerodynamic loads. For example, for a propfan with eight unidentical 
blades with 72 aerodynamic panels for each blade, the total number of panels 
in equations (5) and (6) is 576. 


Unsteady Aerodynamic Model for Frequency Mistuning 


If the vibration mode shapes of all the blades are identical, and if only 
the differences in blade vibration frequencies around the rotor are to be 
accounted for, the computational time and the computer memory required for cal- 
culating the matrices, CAi ] and {G}, can be minimized by calculating the 
aerodynamic loads for each intergroup phase angle independently. For this spe- 
cial case the number of blades in a group, NB, equals one, and the intergroup 
and interblade phase angles are the same. Let the generalized aerodynamic and 
the forcing function matrices for the rth interblade phase angle be [2. r ] 
and {g r } where r = 0,1 ... (NG - 1)). These matrices are calculated by 
using equations (6) and (7) and by paneling only one blade. For yawed flow, 
the matrix {g r } has to be calculated only for r = NG - 1. Then for the 
case of arbitrary frequency mistuning the matrices, [A]] and {G} can be 
expressed in terms of the matrices [& r ]'s and {9NG-1 } following the proce- 
dures in reference 11 as 


where 


CA ] ] 


= CE3 




[ft,] 


C2, NG-1 ] J 


CE3 


-1 


{G} 



{g 0 } 

{g^ 

{g NG-l* 


( 8 ) 


(9) 


CE] = 


CE 0,0 ] [E 0,1 ] • • • 


: NG-1 , NG— 1 


( 10 ) 


■s,r 


i (2irrs/NG) 

= e [I], 


s = 0,1,2 


. NG - 1 


( 11 ) 


9 



Solution of Modal Equations 


For simple harmonic motion 


{q} = {q 0 ) e iut 

Equation (3) leads to the following 


(12) 


- 


1-1 

-</ E- M g g + b K g * - CA 1 3 1 {G} 


(13) 


The amplitude of generalized coordinate vector is calculated for a given 
flow Mach number, yaw angle, rotational speed, and blade setting angle by 
extending the computer program, ASTR0P3, described in reference 5. 


Calculation of Blade Stresses and Strains 

The blade deflected shape and vibration data (modes, frequencies, gener- 
alized masses, and modal stresses) have been obtained by solving equations (1) 
and (2), with NASTRAN code. From the modal stresses and the amplitude of the 
generalized coordinate vector, the elemental stresses have been calculated by 
using the relation 

a 

XXI 
a 

yyi 

a xyi 

where o XX i, ayy^ , and o X yi are the two normal stresses and the shear stress 
of an element of the ith blade in the blade local coordinate system; [S^] 
is the modal stress matrix of an element of the 1th blade; and {q 0 i } is 
the amplitude of the generalized coordinate vector for the ith blade. For 
comparison with stresses measured on the metallic blades, the elemental 
stresses are resolved along and perpendicular to the strain gage direction. 

The resolved stresses of the elements surrounding the strain gages were aver- 
aged to calculate a value at the gage location. For comparison with strains, 
measured on the composite blades, the strains at the gage locations are calcu- 
lated by averaging the resolved elemental strains of the surrounding elements. 
The elemental strains are, in turn, calculated from the elemental stresses and 
the equivalent material properties of the elements. 


- CS i ] «oi> 


(14) 


RESULTS AND DISCUSSION 

The analytical model and the ASTROP3 code were applied to five wind tunnel 
models of 0.62 m (2 ft) diameter which were tested for aeroelastic response in 
yawed flow at NASA Lewis. These models are: (1) SR-5 with ten titanium 

blades; (2) SR-3 with eight titanium blades; (3) SR3C-X2 with eight composite 
blades consisting of 80 percent of the plys along the blade pitch axes and 
20 percent along the +22.5° directions; (4) SR3C-3 with eight composite blades 
consisting of 80 percent of the plys along the blade pitch-axis and 20 percent 
along the +45° directions; and (5) SR3C-X2/SR3C-3 with four SR3C-X2 blades and 
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four SR3C-3 blades mounted alternately (for both frequency and mode shape mis- 
tuning). These five models were experimentally tested for aeroelastic response 
in yawed flow at different Mach numbers, rotational speeds, and blade setting 
angles. The data have been documented in reference 4 for SR-3 and SR-5, and 
in reference 12 for SR3C-X2, SR3C-3 and SR3C-X2/SR3C-3. From this data bank., 
only selected test points were chosen to analyze and validate the analytical 
aeroelastic response model and the ASTR0P3 code. The results are presented in 
eight parts: (1) Analysis cases; (2) calculated and measured nonrotating 
blade frequency table; (3) calculated rotating blade frequency table; (4) cal- 
culated and measured blade stress amplitudes for SR-5; (5) Variation of the 
SR-5 blade stress with the blade stiffness and mistuning; (6) Variation of the 
SR-5 blade stresses with the number of blades (7) calculated and measured 
stress amplitudes for the SR-3; and (8) calculated and measured strain ampli- 
tudes for the SR3C-X2 , SR3C-3 
and SR3C-X2/SR3C-3 . 

In all of the calculations, the frequencies and mode shapes of the blades 
were obtained by using COSMIC NASTRAN. The steady-state deflected position of 
the blade due to centrifugal loading is employed in calculating both vibration 
and aeroelastic response. The structural damping of the blades is not included 
in the calculations. Also for unsteady aerodynamic calculations, each blade is 
discretized with nine (radial) by eight (chordwise) panels. 


Analysis Cases 

Cases considered for analysis are listed in table 1. Four cases each for 
the SR-5 and SR-3 models (covering four rotational speeds, four Mach numbers, 
three blade setting angles and six propfan airflow angles) were chosen. These 
cases cover a wide range of operating conditions. Both the SR-5 and the SR-3 
were chosen for analysis because the blade sweep of the SR-5 model blade is 
larger than that for the SR-3. Also, the number of blades of the SR-5 propfan 
is ten and that of the SR-3 is eight. For the three composite propfan models, 
one test tunnel condition was chosen in order to validate the code with the 
effects of blade ply directions and of both blade frequency and mode shape mis- 
tuning on aeroelastic response due to yawed flow. 


Calculated and Measured Nonrotating Blade Frequencies 

For the five cases considered, calculated and measured nonrotating blade 
frequencies are tabulated in table 2. The average values of the measured fre- 
quencies are listed. Only the frequencies for the first four modes are listed 
in this table since the higher order modes are not important, as will be dis- 
cussed, for forced response due to yawed flow. However, for one of the test 
cases of the SR-5 blade, the effect on blade response of using up to six modes 
is investigated. The measured blade frequencies of each propfan model differ 
slightly from each other because of manufacturing tolerances. Comparing the 
calculated and measured frequencies, the agreement is satisfactory for the 
first two modes and is poor for the third and fourth modes. Although not 
shown, the qualitative agreement between the calculated and measured first and 
second mode shapes is also good. More detailed discussions are presented in 
references 5 and 6. 
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Calculated Rotating Blade Frequencies 


The rotating blade frequencies for the first four modes, for each blade 
and for each test point, are listed in table 3. The frequencies are calculated 
with centrifugal stiffening and softening loads by using COSMIC NASTRAN. The 
effects of hub flexibility are neglected in modeling the blade and, hence, 
there is no structural coupling between the blades. Also, Coriolis forces 
introduced by rotation are not included because they were shown to be negli- 
gible for the propfan blades in reference 13. Detailed discussion on correla- 
tion between calculated and measured frequencies can be found in references 4, 
5, 6, and 8. This comparative table is included in this paper in order to show 
that the excitation frequency due to yawed flow, which is the blade rotational 
frequency, shown in table 3, is placed well below the first mode frequency. 

The importance of this will be discussed in the next section. 


SR-5 Calculated and Measured Blades Stresses 

The calculated stress amplitudes for the four test point cases of the SR-5 
propfan model are compared with the measured ones in table 4. Only the ampli- 
tudes are compared because the interblade phase angle of the calculated stress 
for the tuned rotor is 324° (-36°) which corresponds to a one nodal diameter 
backward traveling wave. The predominant interblade phase angle measured is 
also 324°. Average values of the calculated normal (along the strain gage 
axis) and shear (plane perpendicular to the strain gage axis) stresses of the 
blade finite elements neighboring a strain gage location are tabulated in 
table 4. 

In calculating the blade stresses, a numerical convergence study was con- 
ducted by varying the number of modes for the test point case 8508. This 
study, although the results are not shown here, showed that the major contribu- 
tion to blade stresses comes from the first two modes (about 75 percent from 
the first mode and the remaining from the second mode.) The contribution from 
the higher modes is negligible. This was expected because the excitation fre- 
quency of the forcing function due to yawed flow is the same as the rotational 
frequency and is well below the first mode frequency of the blade at that rota- 
tional speed. Because of this finding, only two normal modes are used for cal- 
culating the stresses listed in table 4. 

The blade stresses are measured at four locations which are shown in fig- 
ure 3. Also included in this figure are strain gage angles measured with 
respect to pitch axis and the radial and chordwise locations of strain gages. 
Gages 1, 5, and 6 measured normal stresses and gage 3 measured shear stresses. 

To give a graphical picture of the correlation between theory and experi- 
ment, calculated and measured stresses for the four test cases are shown in 
figures 4 to 7. The variation of the calculated stress around a strain gage 
is illustrated by a vertical bar in these figures. 

For the test point cases 8508, 8511 and 8610, the correlations in figures 
4, 5, and 7, respectively, are good. The calculated stresses are lower than 
the corresponding measured values. This may be due to the level of blade 
frequency and mode shape mistuning (which usually have adverse effect on 
response), discrepancies in strain measurements, strain gage orientations, or 
the finite element model. The correlation between theory and experiment for 
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the test point case 8607 (fig. 6) is not as good as that in figures 4, 5, and 
7, particularly for strain gage 6. This may be due to the fact that the inflow 
(or yaw) angle is 12° which is significantly higher than 2°, 5°, and 3° for 
the cases 8508, 8511, and 8610, respectively. For this high inflow angle, 
the assumption of attached flow may break down. 


Variation of SR-5 Blade Stresses with Blades Stiffness and Mistuning 

It is well known that the blade stresses are functions of the blade stiff- 
ness and blade mistuning. A decrease in the blade stiffness increases the 
blade stresses. However, the effect of blade mistuning may have either a bene- 
ficial or an adverse effects on blade vibratory response and hence on the blade 
stresses. To illustrate the sensitivity of blade stress due to variation 
in stiffness and type of mistuning, the calculated results of a parametric 
investigation employing the SR-5 propfan test point case 8508 are listed in 
table 5. The model in case 8508-R is the reference rotor with frequencies 
fl = 140 Hz and f£ = 274 Hz, for the first and second modes, respectively. 

For models in cases 8508-1 and 8508-2, the frequencies of all of the blades 
are reduced by 5 percent and are increased by 5 percent, respectively. A 
5 percent reduction in blade frequency increased the blade stress from 1 886x1 0 4 
to 21 96x1 0 4 Nm _ S and a 5 percent increase in blade frequency decreased the 
blade stress from 1886xl0 4 to 1640xl0 4 Nrrr^. This clearly shows that the accu- 
racy of the blade frequency calculation is very important for predicting blade 
stress even at the off-resonance condition. 

The effect of blade alternate frequency mistuning on blade stresses is 
shown in case 8508-3. The frequencies of the odd blades are decreased by 
5 percent and those of the even blades are increased by 5 percent. The blade 
stresses are calculated at the strain gage location 1 of blades 1 and 2. Com- 
paring stresses for the cases 8508-R and 8508-3. a +5 percent frequency mistun- 
ing has increased the blade stress from 1 886x1 0^ to 2181x1 0^ Nnr^ f 0 r blade 1 
and has reduced from 1886x1 0^ to 1651xl0 4 Nm~2 f 0 r blade 2. However, comparing 
stresses for the cases 8508-1 to 8508-3, a +5 percent frequency mistuning has 
not increased the stresses significantly. 

The worst case of mistuning is the one where a "rogue blade" is present. 

In this case, all but one of the blades are identical. To quantify this fact 
for propfan blades, four cases 8508-4 to 8509-7 are considered. The frequen- 
cies of blade 1 are varied from 0.95 to 1.1 times the reference values, fi and 
f2- Comparing the stresses on blade 1 at strain gage location 1, there is a 
significant change in blade stress from 2550xl0 4 Nm - 2 for the case 8505-4, to 
1436x1 O' 4 Nm - 2 for the case 8508-7. This comparison clearly demonstrates the 
importance of quality control in manufacturing the blade and the role of blade 
mistuning in aeroelastic response analysis. 


Variation of SR-5 Blade Stresses with Number of Blades 

Table 6 shows variation of blade stress with number of blades. Evidently, 
the blade stress increases, with a decrease in number of blades. It should be 
noted that with a decrease in number of blades, the interblade phase angle of 
the forcing function is changing even though the frequency is the same. For 
example, the interblade phase angle is 324° for ten blades, 315° for eight 
blades, 300° for six blades and so on. It is also noted from the results that 
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the two lowest damped interblade phase angles for the ten bladed case are 288° 
and 324° which corresponds to two and one nodal diameter backward traveling 
waves, respectively. As the number of blades decreases, the damping of the one 
nodal diameter backward traveling wave mode increases. Also, the amplitude of 
the forcing function increases with a decrease in the number of blades. This 
is because closer blade spacing causes a reduction in blade lift from the 
induced velocity. Thus, the increase in blade response with a decrease in 
number of blades is not surprising, even though it is not obvious. 


SR-3 Calculated and Measured Blade Stresses 

Table 7 compares calculated and measured stress amplitudes for the four 
test point cases of the SR-3 propfan model. For this eight-bladed rotor, the 
interblade phase angle of the calculated stress for the tuned rotor is 315°, 
which corresponds to a one nodal diameter backward traveling wave. The predom- 
inant interblade phase angle measured is also 315°, but there is a slight 
deviation from this angle because of the presence of a small amount of blade 
frequency mistuning. Average values of the calculated normal and shear 
stresses of the blade finite elements neighboring a strain gage location are 
tabulated. 

The blade stresses are measured at five locations which are shown in fig- 
ure 8. Gages 1, 2, 3, and 5 measured normal stress and gage 4 measured shear 
stress. 

For a graphical comparison, the calculated and measured stresses for one 
test case, point 190, are plotted in figure 9. The results show that the cor- 
relation between theory and experiment is good. In the case of SR-5, the cal- 
culated stresses are lower than the measured ones, whereas in the case of SR-3, 
the calculated stresses for some gages are higher and for some gages lower than 
the measured ones. 

Comparing the present correlation between theory and experiment with that 
in reference 3, where a two-dimensional aerodynamic theory was used, the gen- 
eral trend is the same but the present correlation is slightly better in some 
cases, which should be expected from a more refined three-dimensional aerody- 
namic model. But, with a three-dimensional aerodynamic model, the uncertain- 
ties associated with a two-dimensional model, such as the choice of reference 
axis for measuring blade sweep and the validity of sweep correction, do not 
ari se. 


Calculated and Measured Blade Strains for the SR3C-X2, SR3C-3, and 

SR3C-X2/SR3C-3 Propfans 

To validate the present analytical formulation and ASTR0P3 for the predic- 
tion of blade stresses of composite blades with mistuning, three composite 
propfan models were studied. The calculated and measured flutter data of 
these models were presented in reference 6, and the measured aeroelastic 
response data were presented in reference 12. All three models are now ana- 
lyzed for aeroelastic response at one set of wind tunnel conditions: M = 0.36, 

rpm = 5500, Bo.75R * 48°. The calculated and measured strains are compared in 
table 8 at two’strain gage locations which are shown in figure 10. For these 
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cases, strains are compared instead of stresses, because the measured data 
were in this form. 

SR3C-X2 correlation . - The variation of the SR3C-X2 blade natural frequen- 
cies with rotational speed, and a comparison of measured bench natural frequen- 
cies with corresponding calculated ones were presented in reference 5. Since 
the correlation between theory and experiment, as discussed in reference 5, is 
poor for the second mode, three subcases of calculated strains are presented 
in table 8. In the subcase "a", calculated frequencies are used; in the sub- 
case "b", the measured second mode frequency is used; and in subcase "c", the 
''quasi-measured" rotating frequencies are used for the first mode and measured 
nonrotating frequency for the second mode for each blade. The quasi -measured 
frequencies are obtained from the nonrotating measured frequencies with a cor- 
rection to centrifugal stiffening and softening loads. The correction is based 
on Southwell coefficient which is determined from the blade rotating and non- 
rotating calculated frequencies. No correction was applied to the second mode 
frequency because it does not change appreciably with rotational speed. Essen- 
tially, subcase "c" represents a rotor with an arbitrary frequency mistuning. 

Comparing the measured and calculated strains, in subcase "a" for the 
SR3C-X2 model, the measured values at the strain gage location 1 are higher by 
27.8 percent. However, correlation between the measured and calculated strains 
in subcases "b" and "c" is better than subcase "a", with the measured values 
still higher than the calculated by 17 percent. Evidently, this comparison 
shows the significance of the accuracy required in predicting natural frequen- 
cies and the level of frequency mistuning. 

SR3C-3 correlation . - The variation of the SR3C-3 blade measured nonrotat- 
ing frequencies is not available. The value in table 2 represents only one 
blade. Since the correlation between the measured and calculated nonrotating 
first mode frequency (table 2) is not that good, two sets of calculated strains 
are listed in table 8. The subset "d" strain is with calculated first and sec- 
ond mode frequencies, and subset "e" is with quasi-measured first mode rotating 
frequency. As can be seen, the calculated strain for the subsets "d" and "e" 
are greater than the measured value, but the strain for subset "e" is closer to 
the measured value. Again, the calculated results are very sensitive to the 
first mode rotating frequency. 

SR3C-X2/SR3C-3 correlation . - In view of the deviations (table 2) between 
the calculated and measured frequencies for the SR3C-X2 and SR3C-3 blades, two 
sets of calculated results, subcase "f" and subcase "g", are included in 
table 8. These cases represent an alternate mistuned rotor. The strain for 
subcase "f" is based on the calculated frequencies and subcase "g" is based on 
the quasi-measured frequencies. The measured and calculated strains for both 
the SR3C-X2 and SR3C-3 blades in the symmetric group are included. Contrary 
to the results in the subcases "a" and "d" in table 8, the calculated strain 
in subcase "f" is higher than the measured for the SR3C-X2 blade, but is lower 
than measured for the SR3C-3 blade. However, the correlation between theory 
and experiment in subcase "g" with quasi-measured frequencies is better. 

For all of the three models discussed earlier, the correlation between 
calculated and measured strains at strain gage location 2 is very poor. The 
reason for this is not known at this time of writing, and is under 
investigation. 
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For the SR3C-X2 model, the calculated interblade phase angle is 315°, and 
the measured value varied from 312° to 318°; for the SR3C-3 model the calcula- 
ted value is 315°, and the measured value not available; and for the SR3C-X2/ 
SR3C-3 model the calculated value varied from 314° to 316°, and the measured 
value varied from 308° to 317°. The overall agreement between theory and 
experiment is good. 


CONCLUSIONS 

A theoretical aeroelastic response model was developed for predicting 
stresses and strains of metallic and composite propfan blades with blade mis- 
tuning. This model is based on normal modes of the rotating blades and three- 
dimensional lifting surface unsteady aerodynamic theory. Calculated results 
were correlated with published measured strain or stress data for five wind 
tunnel propfan models. The correlations and additional parametric results led 
to the following conclusions: 

1. Calculated and measured stresses or strains agree well in general. 
However, in most of the cases, the calculated values are lower than the meas- 
ured values. Also, calculated and measured interblade phase angle of the 
blade response with and without blade mistuning agree well. 

2. Calculated blade stresses and strains are very sensitive to the first 
mode frequency, as expected, since the excitation frequency is lower than the 
first mode blade frequency in all the cases analyzed. The first two blade nor- 
mal modes are adequate to predict the response due to yawed flow. 

3. Stresses and strains are also sensitive to the blade frequency and the 
mode shape mistuning. A +5 percent alternate frequency mistuning does not have 
a significant adverse effect on blade stresses for the metallic blade (SR-5). 
But both frequency and mode shape mistuning have beneficial effects on compos- 
ite blade strains due to yawed flow. These results demonstrate that both fre- 
quency and mode shape mistuning should be included in aeroelastic response 
analysis. 

4. Calculated results showed that by increasing the number of blades on a 
rotor, a decrease in blade stress amplitudes due to yawed flow occurs. 
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TABLE 1. - ANALYSIS CASE TABLE 


Propfan 

model 

Test point 

Rota- 
ti onal 
speed , 
rpm 

Mach 

number 

Blade setting 
angle, deg 

Tunnel ai r densi ty 

Speed of 

sound 

Infl ow 
angl e, 
deg 





kg m~3 

lb s” 2 ft- 4 

ms“^ 

i n s -1 

SR-5 

8508 

5500 

0.801 

70.8 

0.9393 

1 .8225xl(T 3 

329.8441 

12 986 

2 


8511 

5500 

.797 

70.8 

.9438 

1 .8313 

330.0222 

12 993 

5 


8607 

6000 

.360 

60.8 

1.1883 

2.3057 

339.0900 

13 350 

12 


8610 

6000 

.361 

60.8 

1.1886 

2.3062 

339.0900 

13 350 

3 

SR-3 

190 

8000 

.353 

48.9 

1.2283 

2.3833 

340.4108 

13 402 

8 


273 

8000 

.798 

60.8 

.9010 

1 .9034 

363.5502 

14 313 

2 


277 

7000 

.795 

60.8 

.9820 

1.9053 

363.5502 

14 313 

4 


278 

8000 

.795 

60.8 

.9823 

1 .9060 

363.5502 

14 313 

4 

SR3C-X2 

392 

5500 

.36 

48 

1.1472 

2.2260 

355.3968 

13 992 

8 

SR3C-3 

45 

5500 

.36 

48 

1.1472 

2.2260 

355.3968 

13 992 

8 

SR3C-X2/ 

SR3C-3 

768 

5500 

.36 

48 

1.13744 

2.2070 

350.8248 

13 812 

8 
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TABLE 2. - CALCULATED AND MEASURED NONROTATING BLADE FREQUENCIES 
[Frequencies given in Hz.] 


Propfan 

model 

Mode 1 

Mode 2 

Mode 3 

Mode 4 

Calcu- 

Meas- 

Cal cu- 

Meas- 

Cal cu- 

Meas- 

Cal cu- 

Meas- 


lated 

ured 

lated 

ured 

lated 

ured 

lated 

ured 

SR-5 

mm 

112 

273 

278 


483 

634 

631 

SR-3 


160 

411 

385 


638 

713 

866 

SR3C-X2 

m 

191 

392 

a 371 

700 

693 

752 

987 

SR3C-3 

170 

188 

419 

414 

609 

564 

800 

644 


a Note that this is the measured value before the forced response test. 
Before the earlier flutter test the value is 375 which is reported in 
reference 10. 


TABLE 3. - CALCULATED ROTATING BLADE FREQUENCIES 
[Frequencies given in Hz.] 


Propfan 

model 

Test 

point 

Setting 

angle 

b 3/ 4 R, 

deg 

Rotational 

speed 

Mode 

■ 

2 

3 

■ 

rpm 

Hz 

SR-5 

8508 

70.8 

11 

91.67 

140 

274 

581 

628 


8511 

70.8 


91.67 

140 

274 

581 

628 


8607 

60.8 

6000 

100 

151 

281 

586 

640 


8610 

60.8 

6000 

100 

151 

281 

586 

640 

SR-3 

190 

48.9 

mm 

133.33 

225 

453 

696 

745 


273 

60.8 

Euil 

133.33 

217 

444 

693 

743 


277 

60.8 

Pfll 

116.67 

208 

439 

677 

733 


278 

60.8 

in 

133.33 

217 

444 

693 

743 

SR3C-X2 

392 

48 

5500 

91 .67 

223 

408 

704 

808 

SR3C-3 

45 

48 

5500 

91.67 

215 

435 

654 

822 


TABLE 4. - SR-5 CALCULATED AND MEASURED STRESSES 
[Stresses given in (N/m^)xlO^.] 


Test 

point 

Gage 1 


Gage 3 



Gage 5 


Gage 6 


Measured 

Calculated 

Measured 

Cal culated 

Measured 

Calculated 

Measured 

Cal culated 


Normal 

Normal 

Shear 

Shear 

Normal 

Shear 

Normal 

Normal 

Shear 

Normal 

Normal 

Shear 

8508 

2320 

1886 

295 

201 

168 

43 

183 

311 

248 

999 

383 

243 

8511 

5968 

4701 

737 

301 

418 

109 

687 

775 

618 

2602 

956 

606 

8607 

5142 

3023 

469 

338 

259 

76 

666 

479 

361 

4997 

610 

358 

8610 

1225 

759 

118 

93 

65 

19 

166 

120 

90 

1186 

153 

90 
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TABLE 5. - SR-5 BLADE STRESS VARIATION WITH FREQUENCY MISTUNING 


Case 

Case description 

Blade stress, N/m 2 
( Strai n gage 1 ) 

8508-R 

Test point 8508 

Calculated modes and frequencies (f] and f 2 ) 

1886xl0 4 

8508-1 

Frequencies: 0.95f-j and 0.95f 2 

2196 

8508-2 

Frequencies: 1 . 05f ^ and l.Q5f 2 

1640 

8508-3 

Frequencies: 0.95^ and 0.95f 2 (odd blades) 

l.OSfi and 1.05f 2 (even blades) 

*2181 
b 1 65 1 

8508-4 

Frequencies: 0.9f] and 0.9f 2 (blade 1) 

2550 

8508-5 

Frequencies: 0.95fj and 0.95f 2 (blade 1) 

2167 

8508-6 ! 

Frequencies: 1.05f^ and 1.05f 2 (blade 1) 

1629 

8508-7 

Frequencies: 1-lf] and 1.1 f 2 (blade 1) 

1436 


j*For odd blades. 
“For even blades. 


TABLE 6. - SR-5 BLADE STRESS VARIATION 
WITH NUMBER OF BLADES 


Case 

Case 

description 
(number 
of blades) 

Blade stress, N/m^ 
(Strain gage 1) 

8508-R 

10 

1886xl0 4 

8508-8B 

8 

2041 

8508-6B 

6 

2197 

8508-5B 

5 

2262 

8508-4B 

4 

2305 

8508-1B 

1 

2282 


TABLE 7. - SR-3 CALCULATED AND MEASURED BLADE STRESSES 
[Stresses given in (N/m^JxlO^.] 


Test 

point 

Gage l 

Gage 2 

Gage 3 

Gage 4 

Gage 5 


Calculated 

Measured 

Calculated 

Measured 

Calculated 

Measured 


Measured 



Normal 

Shear 

Normal 

Normal 

Shear 


Normal 

Shear 


Normal 


Normal 



190 

1590 

1273 

323 

732 

623 

55 

1177 

788 

192 

167 

616 

69 

644 

351 

38 

273 

1291 

1553 

354 

594 

699 

66 

956 

914 

205 

136 

680 

60 

523 

372 

31 

277 

2426 

2911 

647 

970 

1264 

103 

1926 

1804 

472 

303 

1339 

154 

1120 

699 

96 

278 

2575 

3075 

700 

1184 

1382 

130 

1905 

1809 

407 

270 

1346 

119 

1042 

736 

61 
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TABLE 8. - SR3C-X2, SR3C-3, AND SR3C-X2/SR3C-3 CALCULATED AND MEASURED BLADE STRAINS 


Test 
poi nt 

Propfan 

model 

Gage 1 

Gage 2 

Measured 

Calculated 

Cal cul ated 
wi th 

measured 

frequency 

Cal cul ated 
wi th 

arbi trary 
mi stuni ng 

Measured 

Calculated 

392 

45 

768 

768 

( SR3C-X2) 
( SR3C-3) 
(SR3C-X2) 
( SR3C-3) 

257x1 0" 6 
195 
143 
177 

<*201xlCr6 

°287 

^199 

f 152 

b 219X10 -6 

e 252 

9164 

9153 

c 2 1 7x1 0“k 

358x1 0 -6 

165 

161 

135 

80x1 O' 6 
b 46 

u? 

'50 


a f, = 223 Hz; f 2 a 408 Hz. 


b f, = 223 Hz; f 2 = 371 Hz. 

c f ! = 223, 230, 219, 219, 223, 223, 225, 230 Hz. 


f 2 = 371, 386, 367, 367, 367, 382, 382, 367 Hz. 



= 215 

Hz; 

f2 

= 435 

Hz. 










= 230 

Hz; 

0 

= 435 

Hz. 










= 223 

Hz; 

f 2 

= 407 

Hz i 

(-X2 

blade); f] 

□ 215 

Hz; 

f 2 

= 435 

Hz 

(-3 blade) 

g fl 

= 223 

Hz; 

f 2 

= 371 

Hz < 

(-X2 

blade); f] 

□ 230 

Hz; 

f 2 

= 435 

Hz 

(-3 blade) 



FIGURE 1. - SR3C-X2 PROPFAN MODEL WIND TUNNEL INSTALLATION. 
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FIGURE 2. - BLADE GROUP SCHEMATIC FOR AN EIGHT-BLADED ROTOR 
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FIGURE 3. - SR-5 STRAIN GAGE LOCATIONS (SUCTION SIDE). 
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FIGURE 4. - COMPARISON OF CALCULATED AND MEA- 
SURED BLADE STRESSES: SR-5 MODEL TEST POINT 
NUMBER 8508. 
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FIGURE 6. - COMPARISON OF CALCULATED AND MEA- 
SURED BLADE STRESSES: SR-5 MODEL TEST POINT 
NUMBER 8607. 
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FIGURE 5. - COMPARISON OF CALCULATED AND MEA- 
SURED BLADE STRESSES: SR-5 MODEL TEST POINT 
NUMBER 8511. 





0 13 5 6 


STRAIN GAGE NUMBER 

FIGURE 7. - COMPARISON OF CALCULATED AND MEA- 
SURED BLADE STRESSES: SR-5 MODEL TEST POINT 
NUMBER 8610. 
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FIGURE 9. - COMPARISON OF CALCULATED AND MEASURED 
BLADE STRESSES: SR-3 MODEL TEST POINT NUMBER 190. 



FIGURE 10. - SR3C-X2 AND SR3C-3 
STRAIN GAGE LOCATIONS (SUCTION 
SIDE). 
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